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Abstract 

We construct holographic superconductors from Einstein-Maxwell-dilaton grav- 
ity in 3+1 dimensions with two adjustable couplings a and the charge q carried by 
the scalar field. For the values of a and q we consider, there is always a critical 
temperature at which a second order phase transition occurs between a hairy black 
hole and the AdS RN black hole in the canonical ensemble, which can be identified 
with the superconducting phase transition of the dual field theory. We calculate 
the electric conductivity of the dual superconductor and find that for the values 
of a and q where a/q \s small the dual superconductor has similar properties to 
the minimal model, while for the values of a and q where a/q is large enough, the 
electric conductivity of the dual superconductor exhibits novel properties at low 
frequencies where it shows a "Drude Peak" in the real part of the conductivity. 
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1 Introduction 



The AdS / CFT correspondence [T] provides an elegant idea to the study of strongly cou- 
pled quantum field theories by relating them to certain classical gravity theories or string 
systems. Nowadays, the holographic correspondence has also become a very efficient 
method to deal with the strongly interacting systems in condensed matter physics. A 
great deal of progress has been made in the application of this holographic method to con- 
densed matter physics. Some nice reviews on this subject can be found in |3l HI O El [7] . 

The high temperature superconductor is an exciting while not completely understood 
subject in condensed matter physics and it remains an unsolved mystery because its 
theoretical basis might be a strongly coupled field theory. It is interesting to use the 
gauge/gravity duality to obtain some insights into the properties of superconductors. 
The simplest model to obtain a holographic superconductor with quite similar behavior 
to real superconductors was first built in [HI |9] through Einstein gravity which is minimally 
coupled to a Maxwell field and a charged complex scalar with a potential term. Below 
some critical temperature Tc, the charged black hole solutions develop a non-trivial hair. 
From the point of view of the dual field theory a U{1) symmetry breaks below at a finite 
charged density because of the condensation of a charged scalar|f| This model naturally 
realized s-wave superconductors and in this paper we will call this model the minimal 
model for simplicity. Complete analysis including the backreactions of this system and 
about the zero temperature limit has been considered in [TH [151 [HI [U] • Later, following 
[l8l [T9l [20] , the fermion spectral function in this system was analyzed and very similar 
behavior to what was seen in the angle-resolved photoemission experiments (ARPES) on 
high Tc cuprates was found in [211 ESI EH]. The realization of p-wave superconductors 
has been studied in [211 [25], and d-wave superconductors in [261 EH [28] . 

Although the minimal holographic superconductor model has achieved much success, 
it is still necessary to consider more generalized non-minimal holographic superconductors 
to incorporate more (or some universal) features of superconductors in real physical 
systems, or even to find the dual description of real superconductors. In [291 EO], a 
general class of superconductors was considered and some universal behavior was found. 
Discussions on other aspects of generalized holographic superconductors can be found 
e.g. in [3ll[32|. 

One of these interesting generalized holographic superconductor models [SU] is the 

Note that, according to the dictionary of AdS/CFT, this U{1) symmetry on the field theory side 
should be a global one, thus the dual field theory is superfluid (TOl [H] ■ Discussions about how to make 
the symmetry a local one can be found in [121 [13]. We assume that this U{1) symmetry will eventually 
be gauged. 
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Einstein-Maxwell-dilaton modeljj For our aim to realize holographic superconductors 
from Einstein-Maxwell-dilaton gravity, there are some additional constraints on the form 
of the action of the gravity theory. One of these constraints is that AdS Reissner- 
Nordstrom (RN) black holes should be solutions to the theory at a finite charge density. 
This excludes the e^'^F'^ type Einstein-Maxwell-dilaton models, as the equation of motion 
for the scalar field 77 in this theory ensures that all the charged black holes carry non- 
trivial hair. This motivates people to consider the Cosh (a?]) type Einstein-Maxwell- 
dilaton models where AdS RN black hole solutions can exist. This was first studied in 
the framework of holographic superconductors in [30] , where it was found that this model 
shares essentially the same physics as the minimal model in [13] when a = 1. In a nice 
early paper |36j, the phase transition between AdS RN black holes and dilatonic black 
holes with neutral dilaton was studied in this type of models and some novel behavior 
in the electric conductivity at low frequencies was found. However, because the dilaton 
considered there is neutral, it does not have a dual superconductor description. In this 
paper, we will consider this model with general values of a and charged dilatons in the 
framework of holographic superconductors. 

For the values of a and charge q we consider, there is always a critical temperature 
at which a second order thermal phase transition occurs between a hairy black hole and 
the AdS RN black hole in the canonical ensemble. Below this temperature the dual 
theory is in a superconducting phase while above this temperature the dual theory is in 
a normal phase. We also study the electric conductivity of the dual superconductor. For 
the values of a and the charge q where a/g is small the dual superconductor has similar 
properties to the minimal model as pointed out in [30]. However, for the values of a and q 
where a/g is large enough, the electric conductivity of the dual superconductor exhibits 
novel properties which are very different from the minimal model at low frequencies, 
e.g. near — )■ 0, a "Drude Peak" arises. This can also be seen from the shape of the 
Schrodinger potential after translating the calculation of the electric conductivity into a 
one dimensional scattering problem. 

In the remainder of this paper, we will first construct the basics of Einstein-Maxwell- 

"'Thc terminology "dilaton" is not so accurate here. Generally "dilaton" refers to a real scalar field 
which is non-minimally coupled to the Einstein-Maxwell theory. Hairy black hole solutions widely 
exist in Einstein-Maxwell-dilaton theory and the holographic dual for Einstein-Maxwell-dilaton model 
is itself very interesting because Einstein-Maxwell-dilaton gravity is very common in the low-energy 
effective theories of string theory and black hole solutions in Einstein-Maxwell-dilaton gravity may 
exhibit some quite special thermal properties |33| . The holography of dilaton black holes has been 
studied in [Ml EH [23 [36] (see also [3ai3i[3i[40l[4Tl[42l[43l[l4]). However, in the framework of 
holographic superconductors the scalar field arising in the bulk theory has to be a complex one as it has 
to be charged under the Maxwell field, and in the title and the remainder of this paper by "Einstein- 
Maxwell-dilaton model" we actually mean the Einstein-Maxwell gravity non-minimally coupled to a 
complex scalar in a similar way to the dilaton field. 
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dilaton model in Sec. 2. In Sec. 3 we give the numerical results of the phase transition 
between hairy black holes and AdS RN black holes in a canonical ensemble. In Sec. 4 
we show the behavior of the electric conductivity of the dual superconductor. Sec. 5 is 
devoted to conclusions and discussions. 



2 Basic Set-up for Einstein-Maxwell-Dilaton Models 

In this section we follow [30] to consider a generalized holographic superconductor model 
built from the most generalized covariant gravity Lagrangian with at most two derivatives 
of fields in 3+1 dimensions. The model has the following field contents: a metric field 
Qp^p, a U{\) gauge field A^, a real scalar field rj and a Stiickelberg field 9, which are 
coupled in the following way 



167rG 



d X\ 



[2.1] 



where G{ri), U{ri) and J{ri) are three functions of the scalar 77, whose forms can affect 
the dynamics of the dual superconductor. 

This system has a f/(l) gauge symmetry and the gauge transformations are the 
standard one — )■ + 9^A, ^ — t- 6* + A, so we can choose the gauge ^ = in 
the following calculations. In fact, we can interpret the scalar field rj as the modulus of 
a complex scalar ip and 9 as its phase, i.e. ip = 776*'^^. 

The Einstein equation of motion for g^y is 

R,. - la^ul^R + l^uiv) - lidavd'-v) - \Giv)F' - Iav)a'^ 

~d,vd.v - \j{^)A,A, - ^G{v)F,,F/ = 0. (2.2) 
The equation of motion for the gauge field is 

vjG{v)F^A-Jiv)A'' = 0, (2.3) 



and the equation of motion for the scalar field 77 is 

The choices of the three functions G{ri), U{ri) and J{ri) are crucial to the building of 
the dual superconductor here. To build a superconductor from AdS/CFT, we need to 
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have an AdS vacuum solution in this system, which requires that U{ri) has a finite and 
positive extremum at r/ = 0. Also we require that AdS RN black hole is a solution to 
this system. As a very simple example, in this paper we consider the following choices 
of the functions G^f]), U{ri) and J{r]) as 



G{r]) = cosh(a?7), 

uiv) = i-\ 

Av) = qW. (2.5) 



2 
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where a, m and q are constants. 

Under this choice, the system has an extra Z2 symmetry: rj — )■ —77. Note that the 
AdS RN black hole would not be a solution to this system if we choose G{ri) to be of 
the form e"'', which can be easily seen from the equations of motion. This choice of 
the three functions (12. 5p has been studied in [30] and [36] for the cases a = 1, g = 3 
and g = with general a, respectively. In [30] it was pointed out that for a = 1 and 
5 = 3, this model gives a dual superconductor which has very similar properties to the 
minimal model studied in [T3]. In [35], the authors found that at g = there are also 
phase transitions for general a and as a increases, some novel properties arise. However, 
as g = in this model, the U{1) symmetry is not broken and the dual field theory does 
not have a superconductor description. In this paper we will analyze this model in detail 
for general values of a and g 7^ and find that as the value of a increases the dual 
superconductor has some different behavior compared to the minimal model discussed 
in [E 



To give the dual superconducting phase, we need a hairy solution with the form 
assumed to be 

9{r) 

A = (t){r)dt, 

T] = Tjir). (2.6) 

The equations of motion can be simplified to be 

X' + + T^.e^Jm' = 0, (2.7) 



2 ' 2g 

^9 

Y2 y' , d,,Gi\ J{r]) 



= 0, (2.9) 



2 G J gGiv) 



V" + V'{--J + '- + ^e>^ W + jTgdvU + ^e>^9,J0^ = 0, (2.10) 
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under the assumption (12. 6p . Note that (12.71) and (12. 8 p are the combinations of the 
Einstein equations of motion for gu and Qrr, while (12. 9p and (I2.10p are the equations of 
motion for and rj, respectively. The equation of motion for g^^ is not independent and 
it can be derived from the four equations abovejf] We can see that the AdS RN black 
hole is a solution to this system with 

X(r) = v{r) = 0, g{r) = ^ " -(^ + f^) + and = p(- - -). (2.11) 

It is difficult to find solutions to the equations of motion analytically and we will do 
this in this paper using numerical methods. We can solve the equations by integrating 
the fields from the horizon r+, which is determined by g{r+) = 0, to infinity numerically. 
There are totally four physical fields which need to be solved: ?7(r), </'(r), x(r) and g{r). 
We demand that (j){r) vanish at the horizon in order for the gauge one-form to be well 
defined at the horizon [8]. At the horizon there are four independent parameters 

r+, V+ = v{r+), E+ = (j)'{r+), x+ = x{r+), (2.12) 

as 5'(r+) = and ?7'(r+) can be determined from the four parameters above using the 
equations of motion expanded near the horizon: 
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r/'(r+)<7'(r+) + ^e^+d,G{v+)El + ^d,U{r,+) = 0. (2.13) 

We can get solutions of the system by integrating the equations of motion given the 
initial values of the four parameters above at the horizon. 

The Hawking temperature for the solution can be calculated as 

r = ^ (l2e-^-/'U{v+) - e^^/'GMElf^ . (2.14) 

Before doing the numerical calculations, we list the three scaling symmetries of this 
system, which can help simplify the calculation. The first one is 

e^^b'^e^, t^bt, (p ^ (p/b, (2.15) 

and we can use this scaling symmetry to set x(^) = at the boundary. The second one 
is 

r br, {t, X, y) (t, x, y)/b, g b^g, 60, (2.16) 



"'Since this is not an obvious observation, we will give a simple proof in the appendix. 
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which can be used to set r+ = 1. The third scahng symmetry is 



r^br, t^bt, £ ^ bi, q ^ q/b, (2.17) 

which rescales the metric to b'^g{r) and A = (j){r)dt to bA. This scahng symmetry can be 
used to set £ = 1 during the calculations. 

At the AdS boundary r — )■ oo the behavior of the fields are the following. For the 
scalar field 

+ (2.18) 

where 

To have a stable theory we need to specify a boundary condition either ip^^'^ = or 
^(A+) ^ -por -9/4 < m'^f < -5/4, either ^^^-^ = or = can be chosen as 

the boundary condition. For —5/4 < m'^i'^, we can only impose the boundary condition 
■ip^^-^ = 0. In this paper we will set m^£^ = —2 for simplicity and A_ = 1 while A+ = 2 
for this value of m. Thus with different choices of boundary conditions we can read off 
the expectation value of a dimension one operator Oi or of a dimension two operator 02- 
For the boundary condition ijj^^^ =0, we can have 

{02)=i^^'\ (2.20) 
and for the boundary condition t/^'-^-' =0, we have 

m = ij^'\ (2.21) 

The boundary behavior of the gauge field is 

0(r) ~ /i - ^, (2.22) 
r 

where fi is the chemical potential of the dual field theory while p is the charge density. 
The boundary behavior of the metric fields x('") and g{r) can be determined from the 
equations of motion to be 

^2 



A^^^r 1 

xvn ~ 

and 



9ir)-^ + ^^^,-— If 1<2A<3, 
2M 

g{r) ~ ^ - — , if 3 < 2A. (2.24) 
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After using the scaling symmetries there are only two parameters at the horizon which 
can be used as initial values: 77+ and E^. At the boundary, we have five parameters which 
give the properties of the dual field theory: /i, p, ip^^\ ip^'^^ and M. Thus by integrating 
out from horizon to infinity, we have a map: 

(r/+,E+)h^(/i,p,^«,^(2)^^)_ (2.25) 



3 Numerical Results for the Condensates 

In this section we give the numerical results for the condensates of the dual superconduc- 
tor in a canonical ensemble, i. e. the charge density of the system is fixed to a value p. For 
the value of m we choose, there can be two boundary conditions for r]{r): = with 
giving the condensate O2 or ip^'^'^ = with ip^'^'^ giving the value of the condensate 
Oi- With the constraints from the boundary condition, the map fl2.25l) reduces to a one 
parameter family of solutions for each choice of a and q. We can think of this parameter 
as being the temperature of the theory at a fixed charge density. 

We find that for the values of a and q we considered, there is always a critical 
temperature below which charged hairy black hole solutions can be found and above 
this critical temperature only AdS RN black hole solutions exist at a fixed nonzero charge 
density. 

The free energy in the canonical system can be calculated from the Euclidean gravity 
action through a Legendre transformation. The free energy of the hairy black hole is 
calculated to be [HI [36] 

i^Hairy = V^(-2M + /ip), (3.1) 

and the free energy of the AdS RN black hole is 

^RN = ^(-rt + ^), (3.2) 

where V is the volume of the (x, y) plane. 

At temperatures lower than the critical temperature Tc, hairy black holes have a 
smaller free energy than the AdS RN black holes. As an illustration we plot the picture 
of free energies for hairy black holes and AdS RN black holes for a = 5 and q = 1 and 
the operator Oi in Figure [TJ Thus below this critical temperature, the system is in a 
superconducting phase while above this critical temperature, the system is in a normal 
phase. 

The origin of this instability of AdS RN black holes at low temperatures is the same 
as in the minimal model [H], which can be attributed to the fact that the effective 
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F 




Figure 1: Free energies for the hairy black hole (the purple line) and the AdS RN black hole 
(the blue line) at a fixe charge density for a = 5,q = 1 and the operator Oi. 



mass of the scalar in the zero temperature limit of the AdS RN black holes violates the 
Breitenlohner-Freedman (BF) bound near the horizon. 

For the AdS RN black hole (l23Til . the temperature is T = (12 - p^)/167r, so the 
extremal limit is at p = 2-\/3 and the near horizon geometry of the AdS RN black hole 
is AdS2X R^, i.e. 

ds"^ = -Q(r -ifdt^ + — - + dx^ + dy'^, = 2V3(r - 1). (3.3) 

6(r — 1)^ 

Plug (13. 3p into (I2.10p and we can recover a wave equation in AdS2 in the r/ ^ 1 limit, 

Tl,ff + - = 0, (3.4) 

where we introduced a new effective mass 

"^eff = ^ , (3-5) 

and a new variable r = r — 1. 

An instability would arise when the mass of t] violates the BF bound in the near 
horizon region, i.e. the AdS2 spacetime, while satisfies the BF bound for four dimensional 
AdS4 spacetime: 

-2g^ -6a^ < > (3.6) 

This kind of instability is very useful to the realization of holographic phase transition, 
and a recent application can be found in |45] . 

The values of the condensates in the superconducting phase as functions of the tem- 
perature for various values of a and q are plotted in Figure [2] for Oi and in Figure [3] for 
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02 respectively. As noticed in j36| [5] , for any given values of a and q there are usually 
several different hairy black hole solutions with the same correct asymptotic behavior 
and we always choose the only one solution with a monotonic scalar profile. 




0.0 0.2 0.4 0.6 0.8 1.0 

Figure 2: Values of the condensate 0\ as a function of the temperature for all the combinations 
of a = 0,1,2,3,5, and q = \ (blue), 3 (black), 5 (green), 8 (purple). For each color, the lines 
from top to down correspond to a = 0, 1, 2, 3, 5 with the same value of q dictated by the color 
respectively. 



In Figure Ej the values of the condensate 0\ as a function of the temperature for 
a = 0,1,2,3,5, and q = 1,3,5,8 are plotted. It can be easily seen that for any fixed 
value of q, the value of Oi decreases as the value of a increases. In Figure |3l the values 
of O2 as a function of the temperature for a = 0, 1, 2, 3, 5, and g = 1, 3, 5, 8 are plotted. 
Note that as the maximum value of O2 for a = and g = 1 is much larger than for the 
other values of a and q, it is not plotted out in Figure El In Figure El it can also be 
discovered that as a increases, the value of O2 decreases for fixed values of q. 



It can be checked that 



VdT 



T=T, 



VdT 



while 



Hairy 



T=T, 



7^ 



T=T, 



VdT^ 



T=T, 



SO the phase transition at is a second order phase transition. We can also see this 



9 



0.0 0.2 0.4 0.6 0.8 1.0 

Wc 

Figure 3: Values of the condensate O2 as a function of the temperature for all the combinations 
of a = 0,1,2,3,5, and q = I (blue), 3 (black), 5 (green), 8 (purple) except a = 0, q = l- For 
each color, the lines from top to down correspond to a = 0, 1, 3, 5, 8 with the same q dictated 
by the color respectively. The case a = 0, q = 1 is not in the figure as it has a much larger 
^/{02)/Tc (about 21) than other cases. 

from the behavior of the condensates near the critical temperature T^. Near Tc, both the 
condensates Oi and O2 behave hke qOi ~ aiT^(l—T /TcY^"^ , for 2 = 1,2 and ai and 02 are 
two constants which differ while a and q change. In Tabled] we list the values of ai and 
02 corresponding to different values of a and q. The behavior qOi ~ ajT^*(l — T/TcY^'^, 
for z = 1, 2 is consistent with the prediction from the mean field theory for second order 
phase transitions [H]. 
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8.0 


6.7 


4.7 


3.4 


2.2 




1 


670 


100 


30 


17 


9.3 


3 


9.0 


8.7 


8.4 


7.5 


5.7 




3 


182 


153 


97 


60 


31 


5 


9.1 


8.9 


8.7 


8.6 


7.6 




5 


155 


146 


126 


98 


57 


8 


9.3 


9.3 


9.3 


9.0 


8.6 




8 


148 


145 


137 


122 


92 



Table 1: Left: The coefficient oi for varfous combinations of a and q. Right: The coefficient 
02 for various combinations of a and q. 

The critical temperature Tc is proportional to ^/p and it also depends on the values 
of a and q. Tc/ y/p increases as a or g increases. In Figure SJ the values of Tc/ y/p as 
functions of a and q are plotted for the operators Oi and O2 respectively. 




Figure 4: Left: The values of Tc/^/p as a function of a and q for the operator Oi; Right: The 
values of Tc/ ^fp as a function of a and q for the operator 



4 Electric Conductivity 

The electric conductivity of the dual superconductor in the superconducting phase can 
be calculated from linear perturbations of gtx and Ar^. around the hairy black hole in 
the gravity side. We consider perturbations with zero momentum: = ax{r)e~^^ and 
gtx = /(r)e~*'^*, and these perturbations can get decoupled from other perturbations. 
The equations of motion for the two perturbations are 

/'--/ + G<P'ax = 0. (4.2) 
r 
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Plug (14. 2p into f l4.ip . and we can obtain a single equation of motion for ax{r): 
'9' X' , d,G A , , f,uj' G(t>'\ , J 



II 



+ <+((^-^y-^K^O. (4,3) 



The asymptotic behavior of the Maxwell field near the boundary is 

a. - af + ^, (4.4) 
and the conductivity of the dual superconductor can be calculated from the formula [H] 



Thus in order to get the electric conductivity, we still have to use numerical calcu- 
lations to get the values of a^^ and a^^ on hairy black hole backgrounds with different 
temperatures for various combinations of a and q. Before performing the numerical cal- 
culations, we need to get the behavior of near the horizon. From (14.31) it can be easily 
seen that should vanish as power law of g{r) near the horizon, i.e. 



-lUJ. 



a^(xg V«^i.=.+ . (4.6) 

The values of electric conductivities as functions of co/T, u/Tc and uj/{q{0-)Y^\ i = 
1, 2 are plotted in Figures El El d El El [lO] for (a, q) = (5, 1), (0, 3), (2, 3), (3, 3), (5, 3), (3, 5) 
respectively. As there is a pole at w = in the imaginary part of the electric conductivity, 
a delta function arises in the real part of the electric conductivity at w = using the 
Kramers-Kronig relations. 

From these figures, we can see that at fixed values of q, the curves of the electric 
conductivities behave similarly to the minimal model (Figure E|) in [E] when a is small, 
and as we increase the value of a some novel behavior of the electric conductivities arises. 
As an illustration, we can compare Figure El with Figure El and find that at the small u 
region, the real parts of the conductivities of the two models behave quite differently. In 
Figure El with a = 5 the real part of the conductivity exhibits a "Drude Peak" in the 
region a; — )■ and has a dip at a finite and small value of u. Among these figures. Figure 
El for a = 5 and q = 1 has the largest deviation from the minimal model and it behaves 
similarly to Figure El qualitatively while it has larger maximum values in the limit a; — )■ 0, 
which are 0.93, 1.82, 7.2, 23 for d and 0.78, 1.44, 3.50, 3.0 for O2 in Figure El This 
novel behavior has also been found in holographic strange metals by [16] and in [36] for 
the case g = 0. 
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Figure 5: a = 5,q = 1. The three pictures on top are electric conductivities for the condensate 
Oi. From left to right, the real part of the electric conductivities are functions of w/T, w/Tc 
and w/q{0'^) respectively. In the middle, the imaginary part of the conductivity is also plotted 
using dashed lines. In the first picture, the lines correspond to to T/Tc = 1,0.801,0.499,0.201 
from top to bottom while in the second and the third pictures, the lines correspond to T/Tc = 
0.499,0.201 from top to bottom. The three pictures on the bottom are for the condensate 02- 
In the first picture, the lines correspond to T/Tc = 1,0.801,0.504,0.200 from top to bottom 
while in the second and the third pictures, the lines correspond to T/Tc = 0.504,0.200 from 
top to bottom. There is a delta function in Re(cj) at w = 0, which is not plotted out. 

As there are peaks for Re((j) at — ?■ for certain values of a and g, the height of the 
peaks depends on the temperature and a, q. In figure [TT] we show the dependence of the 
height on the temperature for some values of a and q. We can see that for a = 5, the 
height does not increase monotonically Uke in the minimal model and has a maximum 
at some finite temperature. This maximum value also decreases as q increases. Thus we 
can see that all these novel properties at small frequencies become more apparent as a/g 
goes larger. 

From Figure im we can also find that as the temperature approaches zero, i.e. T/Tc 
0, the height of the peak vanishes for all the values of a and q considered here. Thus as 
T/Tc 0, the peak vanishes and a gap arises, so the behavior of Re((T) becomes again 
similar to the case of the minimal model. From the Figures El O [71 El [9l [10] we can 
also see that the imaginary parts of the conductivities exhibit poles at w = and have 
minimums at u/Tc ~ 8 for the operator O2, so we can also expect here that the width 
of the gap oJg/Tc ~ 8 as T/Tc 0, which is also the same as in the minimal model 

It is apparent that when we adjust the values of a and q we can get the curves of 
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Figure 6: a = 0,q = 3. The three pictures on top are electric conductivities for the condensate 
Oi. In the first one, the hues correspond to T/Tc = 1,0.792,0.490,0.204 from top to the 
bottom. The three pictures on the bottom are for the condensate 02- The fines in tfie first 
one correspond to T/T^ = 1, 0.804, 0.508, 0.205 from top to bottom. In the four pictures on the 
right, the fines correspond to T/Tc = 0.5,0.2 from top to bottom. 




Figure 7: a = 2,q = 3. The three pictures on top are electric conductivities for the condensate 
Oi. The lines in the first picture correspond to T/Tc = 1,0.78,0.51,0.20 from top to the 
bottom. The three pictures on the bottom are for the condensate 02- The lines in the first one 
correspond to T/Tc = 1,0.84,0.52,0.21 from top to bottom. In the four pictures on the right, 
the lines correspond to T/Tc = 0.5, 0.2 from top to bottom. 
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Figure 8: a = 3,q = 3. The three pictures on top are electric conductivities for the condensate 
Oi. The lines in the first picture correspond to T/Tc = 1,0.8,0.5,0.2 from top to the bottom. 
The three pictures on the bottom are for the condensate 02- The lines in the first picture 
correspond to T/Tc = 1,0.8,0.5,0.2 from top to bottom. In the four pictures on the right, the 
lines correspond to T/Tc = 0.5,0.2 from top to bottom. 




Figure 9: a = 5, g = 3. The top three pictures are for the condensate Oi and the bottom for 
02- The lines in the two pictures on the left correspond to T/Tc = 1,0.8,0.5,0.2 from top to 
bottom. In the four pictures on the right, the lines correspond to T/Tc = 0.5,0.2 from top to 
bottom. 
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Figure 10: a = 3, g = 5. The top three pictures are for the condensate Oi and the bottom are 
for 02- The lines in the two pictures on the left correspond to T/Tc = 1,0.8,0.5,0.2 from top 
to bottom. In the four pictures on the right, the lines correspond to T/Tc = 0.5,0.2 from top 
to bottom. 



2.5 r 




Figure 11: Left: The height of the peak at — > for g = 3, a = (blue), 2 (black), 3 (red), 5 
(cyan) for the operator O2; Right: The height of the peak at cj — )■ for q = 5, = 1 (blue), 3 
(cyan), 5 (red), 8 (black) for the operator ©2. 

electric conductivity with quite different shapes and it is expected that for certain values 
of a and q we can get very close to the shapes of the curves observed in the experiments, 
for example. Fig. 14 in [48] . 

The different small frequency behavior of electric conductivities can also be shown by 
rewriting the equation (12. 9p into a one dimensional Schrodinger equation and expressing 
the electric conductivity using the reflection coefficient. To do this, we introduce a new 
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radial variable u which is defined by 



du 



-dr. 



(4.7) 



At large r, du = dr /r"^ , so u = —1 / r . In this new coordinate system, the boundary is at 
u = and the horizon corresponds to u = — oo. We can define a new field ^ = y/Ga^, 
and fl4.3p becomes a one-dimensional Schrodinger equation in the new coordinate system, 



du"^ 



+ 



UJ 



V{u) 



^ = 0, 



with 



V{u) 



J 



1 d^VG 



(4.8) 



(4.9) 



VG du^ ' 

where the prime denotes the derivative with respect to r. The potential V{u) has an 
extra contribution compared to the case of the minimal model. 

To solve the equation (14.81) with ingoing wave boundary conditions at m = — oo, we 
can first extend the definition of the potential to all u by setting V{u) = for n > and 
then this equation can be solved through a one-dimensional scattering problem. \E'(m) 
can be taken as the wave function and we consider an incoming wave from the right. 
Thus the transmitted wave is purely ingoing at the horizon, which satisfies our boundary 
condition. While at m > 0, the wave function is 



7^e^'^", u>0, (4.10) 



y/Cttx and (14.51) . we can 



(4.11) 



where TZ is the refiection coefficient. Using the definition 
have 

1 -7^ _ I /I dG\ 
"^^"^^ ~ 1 + n '2^\G'd^) 

For the case G = const, the formula above goes back to the original result obtained in 
[T6] . Note that the second term is purely imaginary and will not affect the real part of 
the conductivity. In fact, it vanishes after considering the asymptotical behavior of the 
scalar field r] and we can drop it out from the formula. Thus we can see that the real 
part of the conductivity is fully determined by the refiection coefficient TZ, hence by the 
potential in the Schrodinger equation. 

The shape of V{u) is crucial to the value of the refiection constant and we can analyze 
the property of V{u) near the horizon and the boundary. Using the near boundary m — 
behavior of the fields we can obtain the near boundary behavior of the potential as 



V{u) ~ p^u^ + 



a 



A{2A 



4j 



{Ay 



-u 



|2(A-1) 



(4.12) 
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Thus it vanishes for A > 1, and is a nonzero constant for A = 1 while diverges for 
1/2 < A < 1. Near the horizon u — )■ — oo, the first term in (14. 9p dominates and the 
potential vanishes as V/^e^'^^". Thus in the following, we mainly focus on the behavior of 
V{u) for the operator 02- 

In Figures [12] and [13] we show the shapes of the potential V{u) for some values of 
a and q. In the minimal model, V{u) is always positive along the w-axis and its peak 
becomes wider and higher as the temperature lowers. Thus in the minimal model the 
real part of the conductivity approaches 1 at large frequencies and becomes very small 
at low frequencies. As the temperature lowers, the value of Re((j) at low frequencies also 
gets smaller. In our model, as a/g increases, V{u) can become negative and develop 
a dip on the left of the peak, which is quite different from the shape of V{u) in the 
minimal model. For large frequencies this does not affect the reflection coefficient much, 
thus at large frequencies Re((T) behaves similarly to the minimal model. While at small 
frequencies, the reflection coefficient is affected greatly and thus when a/g is large, Re((T) 
exhibits novel behavior compared to the minimal model. 



V(u)/p V(u)/p 




Figure 12: V(u) for the operator 02- Left: V{u) for q=3 and a = (blue), a = 1 (cyan), 
a = 2 (purple), a = 3 (black), a = 5 (green) respectively, at T/Tc = 0.5. Right: V{u) for a = 3 
and q = I (blue), q = 3 (purple), q = 5 (black), q = 8 (green) respectively, at T/Tc = 0.5. 
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VCu)//o v(u)/p 




Figure 13: V{u) for a = 5,q = 1. Left: for the operator Oi; Right: for the operator 02- Both 
of them are for T/Tc = 1 (blue), 0.9 (cyan), 0.8 (purple), 0.5 (black), 0.2 (green). 



5 Conclusion and Discussion 

In this paper we studied a next-to-minimal realization of holographic superconductors 
from Einstein-Maxwell-dilaton gravity. There are two adjustable constants in our model: 
a and q. a describes the coupling between the dilaton and the Maxwell field and q is the 
charge of the dilaton. 

For the values of a and q we considered, we found that there is always a critical 
temperature at which a second order phase transition occurs between a hairy black hole 
and the AdS RN black hole in the canonical ensemble. Below this temperature the dual 
theory is in a superconducting phase while above this temperature the dual theory is in 
a normal phase. 

We calculated the value of the condensates in the superconducting phase and found 
that as a increases with fixed q, the value of the condensate gets smaller. We also studied 
the electric conductivity of the dual superconductor and found that for the values of a and 
q where a/qis small the dual superconductor has similar properties to the minimal model 
as pointed out in [30]. However, for the values of a and q where a/q is large enough, the 
electric conductivity of the dual superconductor exhibits novel properties which are very 
different from the minimal model at low frequencies: a "Drude Peak" arises at u; — t- and 
the height of the peak depends on the temperature in a non-monotonic way. This can be 
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seen from the shape of the potential V{u) which becomes negative at some finite u. For 
all the values of a and q we considered, the imaginary part of the electric conductivity 
always has a minimum at oj/T^ ~ 8 for the operator O2, which is the same as in the 
minimal model. It is expected that ujg/Tc for the operator O2 is also approximately 8 for 
all the values of a and q. 

One immediate question is about the zero-temperature limit of the superconductor as 
in this paper we mainly focus on the finite temperature behavior. It would be interesting 
to add magnetic fields to this system to study the Meissner effect of the dual super- 
conductor and to understand the structure of fermion spectral functions in this system 
[2T| [22| [23| B9] . It would also be interesting to try to embed these phenomenological 
models into the framework of string theory, just as in [SD], and add DC currents to the 
dual field theory to study the effects of the DC currents as in \5T\ . 
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Appendix 

In this appendix we show that the equation of motion for g^x can be derived from other 
four equations of motion. With our ansatz (12. 6p for the solutions, there are only three 
nonzero equations of motion fl2.2p for g^^, i.e. the equations of motion for g^r, gtt and 
g^x- These three equations of motion are 

- g'^Ru + /^i?.. + 2g^^R,, + |f/(r/) - ^V"r/V.r/ + \G{r,)F^ + \j{r^)A^ = 0(5.1) 
g'^Ru - g'-'Rrr + 2g^^R,, + |f/(r/) + ^V"r/V„r/ + \G{r,)F^ - \j{n)A^ = 0(5.2) 
g'^Ru + r i?„ + ^U{v) - ^V^r/V.r/ - \g{v)F' - i J(r/)A2 = 0(5.3) 

Now our task is to obtain (15.31) from (15.11) and (15. 2p . (15. ip and (15.20 can be simplified to 
be 

g'^Rtt - g'^Rrr = -^V^r^V.r/ + h{r,)A^, (5.4) 
g'^'^Rxx = -lu{ri) - lGiv)F'- (5.5) 



20 



To prove the independence, we need to prove that we can obtain fl5.3p from (15 ■4p and 
For the metric obeying the form of the ansatz (12.61) . we have 



1 + r rx 

9 



+ - 



r{g''Ru - g'^Rr 



(5.6) 



Substituting (15. 4p and (15.51) into the above equahty (15. 6p . and using (12. 4p and the 
following useful expression 



2 x' 

drd^f] = V Tj — i )d^'n) 

r 2 



we can obtain 



-lG{v){Fy + \j{ri){A'r + jjJ{v)A'{g' - gx') 



(5.7) 



(5.8) 



From the equation of motion for the gauge field (12.31) . i.e. 

(-y + ^)G(r/)F^* + d^G{r])7]'F'-' + G{r])drF^' = JA\ 

we obtain the following identity on shell 

--G{r])F^ - 2d,G{r]WF' + U{r])FrtA\ 
r 



GmF^y- 

and at the same time, we also have the following off- shell relation 



{A^y = 2FrtA' + ix' - ^-)A'. 

9 



(5.9) 



(5.10) 



(5.11) 



Substituting (I5.10p and (15. lip into (15. 8p . we obtain (15. 3p . Thus only two equations of 
motion in (12. 2p are independent. 
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